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1 Introduction 


Ce document présente des corrections en Xcas du stage d’algorithmique pro- 
posé en 2009/10 par l’IREM de Grenoble. Il est distinct du document Solutions 
de l’IREM de Grenoble, dont les auteurs n’ont pas souhaité intégrer de corrections 
programmées en Xcas. Nous profitons de l’indépendance de ce document pour 
proposer dans certains cas des algorithmes et programmes que nous jugeons mieux 
adaptés (utilisant d’autres structures de données, ou des fonctions avec paramètres 
et valeur de retour au lieu d’avoir à faire des entrées et des sorties fastidieuses, 
ou en simplifiant certains algorithmes) et pour faire le lien avec des algorithmes 
classiques. Les thèmes abordés sont les suivants 

— l'affectation, 

— les listes et les chaînes de caractères avec Xcas, 

— le maximum et le minimum d’une liste de nombres, 

— la factorielle 

— les algorithmes de PGCD, dont l’algorithme d’Euclide. Pour l’algorithme 
(inefficace) utilisant la décomposition en facteurs premiers, on présente le 
crible d’Eratosthene, 

— le crible d’Eratosthene, 

— les changements de base, où nous utilisons des chaînes de caractères ou des 
listes de nombres pour écrire les entiers dans une base différente de 10. On 
y fait le lien avec l’algorithme de Hôrner d’évaluation d’un polynôme. 

— la somme de grands entiers, le produit par un chiffre puis le produit de deux 
grands entiers en utilisant les fonctions précédentes, et en évoquant le lien 
avec les polynômes. 

— le jeu de recherche d’un nombre par une méthode dichotomique, on y fait le 
lien avec la recherche de 0 d’une fonction continue 

— les algorithmes avec sortie graphique sont abordés avec la tortue de Xcas 
mais aussi avec les fonctions géométriques de Xcas 

On pourra aussi consulter avec profit des documents complémentaires dont 

— Le tutoriel algorithmique de programmation Xcas paru dans sesamath, novem- 
bre 2009 

— Le manuel de programmation de Xcas (dans le menu Aide, Manuels) 

— Le tutoriel Débuter en calcul formel (menu Aïde) 

— les documents de la page pédagogique de Xcas 
http://www-fourier.ujf-grenocble.fr/-parisse/algo.html 


2 Affectation 


2.1 Échange du contenu de 2 mémoires 


Pour faire un échange du contenu de deux variables on peut écrire la fonction 
echange qui utilise une variable locale c qui sert à stocker la valeur de à, puis on 
met dans a, la valeur de b et dans b, la valeur de c. On tape : 


echange(a,b) :={ 
local c; 
C:=a; 
a:=b; 


b:=c; 
retourne a,b; 
hi; 
Puis on tape 
a,b :=echangel(a,b) 
Remarque : On peut aussi taper directement : 
a,b :=b,a 


Cela n’est pas identique à la suite d’instructions : a :=b;b :=a; 
Exercice Comparez : 

a :=5;b :=3; 

a :=atb;b :-=a-b; 

on obtient alors 8 dans a et avec 5 dans b 

avec : 

a :=3;b :=5; 

a,b :=a+tb,a-b; 

on obtient alors 8 dans a et avec 2 dans b 
Exercice Echanger les contenus de a et b sans utiliser de variable intermédi- 

aire. Donner une preuve de la méthode utilisée. 


Solution : 
a:=a+b; 
b:=a-b; 
a:=a-b; 


La preuve faite par le calcul formel. On tape dans un éditeur de programmes : 


a:=X; 

b:=y; 
afficher(a,b); 
a:=atb:; 
b:=normal(a-b); 
a:=normal(a-b); 
afficher(a,b); 


On valide avec | OK | ou avec la touche F9 et on obtient : en bleu les 2 affichages 
intermédiaires : x, y 


YrX 
et comme réponse 


x,Y/1,Done,x,y,1 


où (1 est le résultat d’un affichage et Done le résultat de a :=a+b : ; (car cette 

instruction se termine par : ;). 

Remarque : cet algorithme n’a aucun intérêt en pratique, on a vu comment échanger 
le contenu de 2 variables en une ligne avec le langage de haut niveau de Xcas. Avec 

des langages de bas niveau, il existe en général une instruction d’échange (par ex- 

emple swap en C++), et c’est aussi le cas en assembleur (par exemple XCHG sur 

architecture x86). 


2.2 Permutation circulaire de 3 mémoires 


On tape dans un éditeur de programmes : 


a: —X}7 

b:=y; 

C:—Z7y 
afficher(a,b,c); 
d:=a:; 

a:=b; 

b:=c; 

c:=d; 


afficher(a,b,c); 


On valide avec | OK | ou avec la touche F9 et on obtient : en bleu les 2 résultats des 
affichages : x, v,z et 

V, Z, x et Comme réponse 

X/YrZrl1,Done,v,2z,x,1 (1 est le résultat des affichages et Done le résultat 
de d :=a : ; car cette instruction se termine par : ;). 

ou bien on tape : 


a:=X; 
b:=y; 
C:=Z; 


afficher(a,b,c); 
a,b,c:=b,c,a; 
afficher(a,b,c); 


On tape dans un éditeur de programmes : 


a: —X7y 

b:=y; 

CZ; 
afficher(a,b,c); 
d:=c:; 

c:=b; 

b:=a; 

a:=d; 


afficher(a,b,c); 


On valide avec | OK | ou avec la touche F9 et on obtient : en bleu les 2 affichages 
intermédiaires : 

XyYr2Z 

Zr XI YŸ 

et comme réponse 


X,YVrZr1,Done,2z,x,y,1 


où 1 est le résultat des affichages et Done le résultat de d :=a : ; (car cette 
instruction se termine par : ;). On peut permuter les 3 variables en une seule 
instruction a,b,c :=c,a,b; qu’on peut prouver par : 


a: —X}y 

b:=y; 

C:—Z7y 
afficher(a,b,c); 
a,b,c:=c,a,b; 
afficher(a,b,c); 


3 Les listes et les chaînes de caractères avec Xcas 


3.1 Les listes 


Une liste de nombres est une suite de nombres séparés par une virgule que l’on 

met entre des crochets. 
Exemple 

— []} est la liste vide (dim ( [ ] } =0), 

— L :=[0,1,10,16] est une liste de 4 éléments (dim (L) —4) 

— le premier élément est L[0] et il vaut 0... 

— le dernier élément est L [3] et il vaut 16. 

— Siontape:L[2] :=4; alors Lestla liste [0,1,4,16]. 

— [2 $4] désigne la la liste [2,2,2,21],[k$ (k=0..3) ] désigne la liste 
[0,1,2,3]et [k”2$(k=0..4)] désigne la liste [0,1,4,9,16], 

— si festune fonction, makelist (f,0,3) désignelaliste [F(0),£(1),£f(2),£(3)] 
parexemplemakelist (x->x, 0, 3) désignela liste [0,1,2,3]etmakelist(x->x"2,0,4) 
désigne la liste [0,1,4,9,16]. 

On peut taper une fonction Permutation qui effectue une permutation sur les 
éléments d’une liste qui consiste à mettre le premier élément en fin de liste (tail (L) 
est la liste L privée de son premier élément) et append (L, a met a à la fin de la 
liste L alors que prepend (L, a met a au début de la liste L. On tape : 


Permutation(L):={ 
retourne append(tail(L),L[O]l); 
}:; 


3.2 Les chaînes de caractères 


Une chaîne de caractères est délimitée par deux ". 

Exemple 
— "est une chaîne vide (dim ("")=0) 
— $S :="salut" est une chaîne de 5 caractères (dim (S)=5) 
— le premier élément est S [0] et il vaut le caractère "s".… 
— le dernier élément est S [4] et il vaut le caractère "t". 
— C :="256" est une chaîne de 3 caractères (dim (C) =3) 
— le premier élément est C [0] et il vaut le caractère "2"... 
— le dernier élément est C [2] et il vaut le caractère "6". 
— siontape:C[2] :="4";,alors C est la chaîne "254". 
Pour concaténer 2 chaînes on peut utiliser +. 

Ontape:B :="Bon"+"jour" 

et Best la chaîne "Bonjour". 


Lorsqu'on a une chaîne de chiffres on peut la convertir en un nombre avec 
expr. On tape : expr ("136"). On obtient le nombre 136. 
Réciproquement, on peut convertir un nombre en une chaîne de chiffres soit avec 
string, soit avec la concaténation avec la chaîne vide. 
Ontape:a :=136;b :=string(a) ;c :=a+"";d :=""+a 
On obtient la chaîne "136" dans b, dans c et dans qd. 


4 Maximum et minimum d’une liste de nombres 


4.1 Le maximum d’une liste de nombres 
On tape : 


//renvoie le plus grand nombre de L 
Max (L) :={ 
local M, j,d; 
d:=dim(L)-1; 
si d==0 alors 
retourne "liste vide" 
fsi; 
M:=L[0]; 
pour j de 1 jusque d faire 
si M< Lf]j] alors 
M:=L[j]; 
EST; 
fpour; 
retourne M; 


}:; 


4.2 Le maximum et le minimum d’une liste de nombres 
On tape : 


//renvoie le plus grand et le plus petit nombre de L 
Maxmin(L):={ 
local M,m, j,d; 
d:=dim(L)-1; 
si d==0 alors 
retourne "liste vide" 


fsi; 
M:=L[0]; 
m:=L{[0]; 


pour j de 1 jusque d faire 
si M< Lf]j] alors 
M:=L[j]; 
sinon 
si m> L[]j] alors 
m:=L[j]; 
fsi; 


EST; 
fpour; 
retourne M,m; 


}:; 


4.3 Les deux plus grands nombres d’une liste de nombres 
On tape : 


// renvoie les 2 plus grands nombres de L 
// par ordre decroissant 
Max2 (L):={ 
local M1,M2, j,d; 
d:=dim(L)-1; 
si d==0 alors retourne "liste vide" fsi; 
si d==1 alors retourne L[0],L[0];fsi; 
Mi:=max(L[O0O],L{[1]); 
M2:=min(L{[0],L[1]); 
pour j de 2 jusque d faire 
si Mi1< L[j] alors M2:-=M1;M1:=Lfjl; 
sinon 
si M2<L[]j] alors M2:=L[j]; fsi; 
fsi; 
fpour; 
retourne M1,M2; 


} 


5 Factorielle 


5.1 Factorielle et itération 


On tape : 
//fonction qui renvoie n! si n>=0 
Factorielle(n) :={ 

local j,fact; 

fact:=1l; 


pour j de 1 jusque n faire 
fact:-=factx;; 
fpour; 
retourne fact; 
}:; 


5.2 Factorielle et récursivité 


On tape : 


Factorielleri(n) :={ 
si n==0 alors retourne 1l;fsi; 


retourne nxFactorieller(n-1); 


Le 


| 1 
Exercice Calculer 1 - T + DT 
On tape : 


Sommee (n) :={ 
local fact, j; 
fact:=1; 
S:=1; 
pour j de 1 jusque n faire 
fact:-=factx;; 
S:=S+1/fact; 
fpour; 
retourne evalf(S),s:; 
}:; 


On calcule ici dans la foulée la valeur de j! (on n’utilise pas la fonction précédente 
pour ne pas ralentir l’exécution du programme). On tape : 
evalf(e),Sommee(20) 

On obtient : 

2.71828182846, 
2.71828182846,6613313319248080001/2432902008176640000 
On remarque que Xcas fait du calcul exact et que S est une fraction alors que 
evalf(S) renvoie une valeur approchée de S&. 

Si vous changez le nombre de Digits par exemple : 

On tape : 

Digits :=20 : ;evalf(e),Sommee(20) 

On obtient : 
"Done",2.71828182845904523535,2.71828182845904523532, 
6613313319248080001/2432902008176640000 


6 Pgocd 


6.1 Pgcdetitération 


L’algorithme d’Euclide se sert des deux propriétés : 
— sib==0 alors pgcd(a,0)=a 
— sib !=0 alors pgcd(a,b)=pgcd(b, r) sia=bxq+r avec 0O<=r<b 
Dans Xcas on a q=iquo(a,b)etr=iremi(a,b). 
Donc le pgcd de a et b est le dernier reste non nul dans les divisions successives de 
a par b, puis de b par r etc... 
Voici le programme de l’algorithme d’Euclide en utilisant la commande irem(a,b) 
qui renvoie le reste de la division euclidienne de l’entier a par l’entier b. 


//algo d'Euclide 
//a et b sont des entiers positifs ou nuls 
pgcd(a,b):=1{ 

local r; 


tantque b!=0 faire 
r:=irem(a,b}; 


a:=b; 
bi=Tr; 
ftantque; 


retourne a; 
}:; 
Si a est plus petit que b, la première itération de la boucle échange a et b, il n’est 


donc pas nécessaire de faire un test avant la boucle. 
Sans utiliser la commande irem(a,b), on utilise l’instruction : 


tantque a>=b faire a :-=a-b ftantque 


qui met dans a le reste de la division euclidienne de l’entier a par l’entier b. On 
tape : 


//le pgcd avec l’algo d’Euclide sans utiliser irem 
pgcde (a,b) :={ 
tantque b!=0 faire 
tantque a>=b faire a:-=a-b; ftantque; 
a,/b:=b,a; 
ftantque; 
retourne a; 
}:; 
On peut aussi écrire : 


//le pgcd avec l'algorithme des differences 
pgcdd(a,b):=1{ 
si b==0 alors retourne a fsi; 
si a==0 alors retourne b fsi; 
repeter 
tantque a>b faire a:-=a-b;ftantque; 
tantque b>a faire b:=-b-a;ftantque; 
jusqua b==a; 
retourne a; 


ir 


6.2 Pgcd et récursivité 


On a: 
pgcd(a,0)=aet 
pgcd(a,b)=pgcd(b,r) sia=bxq+r avec O<=r<b 
Dans Xcas onar=irem(a,b) etq=iquo(a,b). 
On tape : 


//le pgcd avec l'algorithme recursif 
pgcdr(a,b) :={ 
Si b==0 alors retourne a; fsi; 
retourne pgcdr(b,irem(a,b)); 
}:; 


10 


On tape : 
pgcd(32,424),pgcde(32,424),pgcdd(32,424),pgcdr(32,424) 
On obtient : 

8,8,8,8 


7 Le crible d’Eratosthene 


7.1 Description 


Pour trouver les nombres premiers inférieurs ou égaux à N : 
1. On écrit les nombres de 2 à N dans une liste. 
2. On met 2 dans la case P. 


3. SP X P < N il faut traiter les éléments de P à N : on barre tous les 
multiples de P à partir de P x P. 


4. On augmente P de 1 
Si P x P est strictement supérieur à N, on arrête 


5. On met le plus petit élément non barré de la liste dans la case P. On reprend 
à l’étape 3 


7.2 Écriture de l'algorithme avec Xcas 


//renvoie la liste des nombres premiers<=n selon eratosthene 
criblei(n):={ 

local tab,prem,p, j; 

//on ecrit dans trab la liste des entiers de Oan 

tab:=-makelist(x->x,0,n); 

//1 n'est pas premier donc on met 0 ds tabl1] 

tab[1]:=0; 

p:=2; 

//on barre les multiples de p qui sont compris entre pxp et n 

//en mettant 0 dans la case correspondante 

tantque (pxp<=n) faire 

pour j de pxp jusque n pas p faire 


tab[j]:=0; 
fpour; 
p:=-p+l; 
//on cherche le premier element non barre dans tab 
tantque ((pxp<=n) et (tab[p]l==0)) faire 
p:=p+l; 
ftantque; 
ftantque; 
//on remplit la liste prem avec les elements non nuls de tab 
prem:={[]; 
pour j de 2 jusque n faire 
si (tab[]j]l!=0) alors 
prem:=append(prem, j); 
fsi 
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fpour 
retourne (prem) ; 


}:; 


On tape : 

crible(64) 

On obtient : 
[2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61] 


8 Changement de base 


8.1 Écriture en base b d’un nombre entier avec 2<=b<=10 à l’aide 
d’une liste 


On écrit un nombre n en base b avec des “chiffres” bo, …, b.. Ces chiffres seront 
représentés par des entiers tels que 1 < bo < b,0 < by < b,...,0 < b, < b: 


n = bob° + bb + ...b,_1b + b, 


On les stockera dans une liste de nombres entiers B. Avec Xcas on a pour j = 0..s, 
b; —BT;]. 

On remarque que : 
irem(n,b)=b, et 
iquo(n,b)=Qg = bob°=1 + b1b5—2 + ...b,_1 
D'où irem(q,b)=b,_1 etc … 

On va donc faire une boucle qui calculera irem(n,b) puis qui recommencera 
avecn :=iquo(n,b). 

Il faut faire attention pour générer la liste car on obtient en premier son dernier 
élement, donc il faudra mettre à chaque êtape l’éément suivant devant ceux que l’on 
a déjà obtenu.….mais heureusement dans Xcas, il y a la commande prepend(L, a) 
qui renvoie une liste où lélément a est mis au début de la liste L. 

On tape : 


//ecriture en base b avec une liste B d’entiers de 0..b-1 
dix2B(n,b) :={ 
local B; 
si n==0 alors retourne [0] fsi; 
B:=[]; 
tantque n>0 faire 
B:=prepend(B,iremi(n,b)); 
n:=iquo(n,b); 
ftantque; 
retourne B; 


Fi 


On tape : 
dix2B(139,7) 
On obtient : 
[2,5,6] 
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et on a bien : 6 +5 x 7 + 2 x 49 — 139 

Remarque 

On tape : 

dix2B(139,10) 

On obtient : 

C1,:3,91 

On remarque que dix2B peut servir à transformer un nombre entier en la liste des 
nombres servant à son écriture. 


8.2 Écriture d’un nombre entier lorsqu’on donne la liste de ses coef- 
ficients en base b avec 2<=b<=10 


On entre B = [bo, b1...b,] avec pour j = 0..s 0 <= b; < b et on veut trouver 
n tel que : n = bob° + b1b51 + ...b,_1b + b.. 
Remarque 
n est la valeur du polynôme P(x) = box + bi! + ...b, 1x + b, en x = b. 
On remarque, par exemple pour s = 4 que : 
n = bob + b1b% + bb? + bad + ba = ((((0 x b + bo)b + b1)b + b2)d + b3)d + ba 
On initialise n à 0 et on calcule n en lui ajoutant B [0], puis en le multipliant par 
b, puis, on lui ajoute B[1] etc. 
On tape : 


//n est ecrit en base b avec une liste 
//pour calculer n en base 10 on retrouve l’algo de Horner 
B2dix(B,b) :={ 

local s,j,n; 

s:=dim(B)-1; 

n:=0; 

pour j de O0 jusque s faire 

//si B[j]>=b alors retourne "erreur" fsi; 


n:=nxb; 
n:=n+B[;j]l; 
fpour; 


retourne n; 


}:; 


On tape : 

B2dix([2,5,6],7) 

On obtient : 

139 

etona bien: 6+5%x7+2 x 49 = 139 
Remarque 

On tape : 

B2dix([1,3,9],10) 

On obtient : 

139 

B2dix([1,3,29],10) 

On obtient : 

159 
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On remarque que B2dix peut servir à transformer une liste de nombres considérés 
commes les coefficients d’un polynôme en un nombre entier qui est la valeur du 
polynôme pour x = le 2-ième argument. C’est pourquoi nous commentons dans 
B2dix le test : 

//si B[j]>=b alors retourne "erreur" fsi; 

Pour savoir si la liste B repésente bien une écriture en base b, on peut écrire : 


estbienecritl(B,b}):=f{ 
local j,s; 
s:=dim(B)-1; 
pour j de O0 jusque s faire 
si B[j] >= b alors retourne 0; fsi; 
fpour; 
retourne l; 


}:; 


8.3 Écriture en base b d’un nombre entier avec 2<=b<=10 à l’aide 
d’une chaîne de caractères 


On tape : 


// ecriture de n en base b dans une chaïine 
dix2bci(n,b):={ 
local C; 
Eee 
tantque n>0 faire 
C:=irem(n,b)+C; 
n:=iquo(n,b); 
ftantque; 
retourne C; 
}:; 
On tape : 
dix2bc(139,7) 
On obtient : 
V2:5:6" 


84 Écriture d’un nombre entier lorsqu’on donne ses chiffres en base 
b avec 2<=b<=10 sous la forme d’une chaîne de caractères 


Pour savoir si la chaîne B repésente bien une écriture en base b, on peut écrire : 


estbhbienecritc(B,b):=f{ 
local j,s; 
s:=dim(B)-1; 
pour j de O0 jusque s faire 
si expr(B[j]) >= b alors retourne 0; fsi; 
fpour; 
retourne l; 


}:; 
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Pour convertir, on tape : 


//nc est ecrit en base b avec une chaine 
bc2dix(nc,b):={ 
local s,j,n; 
s:=dim(nc)-1; 
n:=0; 
pour j de O0 jusque s faire 
// si exprinc[j])>=b alors retourne "erreur" fsi; 


n:=nxb; 
n:=n+texpr(nc{jl); 
fpour; 


retourne n; 


F2 


On tape : 

bc2dix( "256",7) 
On obtient : 

139 


8.5 Nombres écrits sans 2 en base 3. 
8.5.1 L’énoncé 


On veut afficher en base dix la suite ordonnée des entiers dont l’écriture en 
base trois ne comporte que des 0 ou des 1 (pas de 2). 


1. Calculer à la main les huit premiers termes de cette suite. 
2. Décrire un algorithme qui donne les 128 premiers termes de cette suite. 


3. Écrire une fonction qui renvoie la liste des n premiers termes de cette suite. 


8.5.2 La correction 


1. Voici les 8 premiers termes de cette suite : 

[0,1,3,4,9,10,12,13] dont lécriture en base 3 est : 

CCOT SLT TROT: LL T1 00:07 071,2C1%0, 11; 613-101. [1;1»11] 

2. Il y a plusieurs algorithmes possibles : 

— On écrit les nombres de 0 à N en base 3 et on met dans la liste réponse 
ceux qui ne contiennent pas de 2 dans leur écriture en base 3. On s’arrete 
quand la liste réponse contient n éléments, 

— On écrit les nombres de 0 à n en base 2 et on considère cette écriture 
comme étant celle d’un nombre écrit en base 3. 

— On regarde comment les termes de la suite se déduisent les uns des autres. 
On peut remarquer que : 
0=3*0, 1=3*0+1, 3=3*1, 4=3*1+1, 
9=3*3, 10=3*3+1, 12=3*4, 13=3*4+1... 

— On regarde comment sont faits les termes de la suite. On peut remarquer 
que 
— 3=0+3, 4=1+3 donc 3 et 4 ont èté obtenu à partir de O et 1 en leur 

ajoutant 3, 


15 


— 9=0+9, 10=1+9, 12=3+9, 13—4+9 donc 9,10,11 et 12 ont èté obtenu à 
partir de 0,1,3 et 4 en leur ajoutant 9... 

Les 8 prochains termes de cette suite seront obtenus en ajoutant 33 à cha- 

cun des 8 premiers. On obtient ainsi les 16 premiers termes de cette suite. 

Puis les 16 prochains termes de cette suite seront obtenus en ajoutant 34 à 

chacun des 16 premiers termes etc... 


3. On va traduire avec Xcas le premier et le dernier algorithme décrit ci-dessus. 

Voici les fonctions de Xcas que l’on va utiliser : 

— La fonction est_element (L,a) qui teste si a est dans la liste L et qui 
renvoie O0 ou n+1 si nest l’indice de la première occurence de à dans L, 

— La fonction convert (n,base,b) (respconvert (L,base,b))qui 
convertit un entier n en la liste des coefficients en base b dans l’ordre 
croissant (resp qui convertit la liste L des coefficients en base b dans l’odre 
croissant en un entier n). On peut aussi utiliser les fonctions dix2B (n, 2) 
etB2dix(L,3) définies précédement mais les calculs seront plus lents. 

— Pour ajouter un nombre à à chacun des termes d’une liste L on peut 
utiliser l’addition de L et de la liste [a,a..a] de la taille de L et écrire 
L+[a$dim(L)] 

Voici les programmes correspondants aux algorithmes 1 et 4 décrits précéde- 

ment à la question 2. (Pour plus de détails voir le manuel d’algorithmique de 

Xcas). 

On tape les fonctions pasde21 (n) et pasde24 (n) qui renvoient les n 

premiers termes de la liste demandée. La variable p contient à chaque étape 

la dimension de la liste I. 


pasde21i(n):={ 
local L,}),J,p; 


L:=[0]; 
p:=1; 
j:=1; 


tantque p<n faire 
J:=convert (j,base,3); 
//J:=dix2B(3,3); 
si not(est_element (2,J)) alors 
L:= appendi(L, j); 
p:=p+l; 
fsi; 
j:=3j+1; 
ftantque; 
retourne L; 


} 


— La variable j contient le nombre d’iterations : à chaque étape on a p — 27 
et puis3j= 3. À la fin de la boucle tantque, L a 2 — p>=Nn 
éléments : il faut donc raccourcir la liste L(L :=mid(L,0,n)) 
pasde24(n):={ 
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local L,j,p,puis3;; 


L:=[0]; 
j:=0; 
p:=1; 
puis3j:=1; 


tantque p<n faire 
L:=concat (L,L+[puis3j$pl); 
//L:=concat (L,map(L,x->x+puis3))); 


j:=3+1; 
puis3j:-3+puis3); 
p:=2*xp; 

ftantque; 


L:=mid(L,0,n); 
retourne L; 


} 


Dans le programme ci-dessus, on calcule des termes pour rien... On mod- 
ifie donc le programme : 
pasde24bi{(n):={ 

local L,]j,p,puis3j; 


L:=[0]; 
j:=0; 
p:=1l; 
puis3j:=1; 


tantque 2xp<=n faire 
L:=concat (L,L+[puis3j$pl); 
Jr; 
puis3j:-3+puis3;; 
p:=2*xp; 
ftantque; 
L:=concat (L,mid(L,0,n-p)+[puis3j$(n-p)l);; 
retourne L; 
}:; 
Le dernier algorithme est le meilleur... 


9 Somme, Produit 


9.1 Nombre de chiffres d’un entier 
On tape pour ne pas utiliser la fonction 1n : 


taillei(n):={ 

local t; 

si n==0 alors retourne 1; fsi; 

t:=0; 

tantque n!=0 faire 
n:=iquo(n,10); 
t:=Et+1; 

ftantque; 
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retourne t; 
}:; 


ou encore on peut se servir de dix2B et on écrit : 
taille(n) :-=dim(dix2B(n,10)) 

On tape : 

taille(123456789) 

On obtient : 
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9.2 Somme de 2 grands entiers 


On va écrire la fonction somme de deux entiers a et b comme on le fait à la 
main. Pour cela, on va prendre comme paramètre de la fonction Somme les listes À 
et B constituées par les nombres servant à l’écriture de ces deux entiers (on pourra 
transformer chacun des 2 nombres a et b en ces listes avec dix2B et obtenir À et 
B). 

On complète la plus petite liste en mettant des Zéros au début de façon à avoir 2 
listes de même longueur. 

Au début il n’y a pas de retenue donc ret :=0. On remplit la liste Result, qui 
contiendra le résultat final, en commençant par le chiffre des unités c’est à dire par 
la fin de la liste. 

On tape : 


//on entre des listes de nombres compris entre 0 et 9 
//on renvoie une liste de nombres compris entre 0 et 9 
Somme (A,B) :={ 
local j,sa,sb,s,c,ret,Result; 
Result:=[]; 
sa:=dim(A)-1; 
sb:=dim(B)-1; 
si sa>sb alors 
B:=concat ([0$ (sa-sb)],B);s:=sa; 
sinon 
A:=concat ([0$ (sb-sa)],A);s:=sb; 
fsi; 
ret:=0; 
pour j de s jusque 0 pas -1 faire 
c:=A[j]+B[jl+tret; 
Result[jl:-=irem(c,10); 
ret:=iquo(c,10); 
fpour; 
si ret>0 alors Result:=-prepend(Result,ret); fsi; 
retourne Result; 


On tape par exemple : 


A :=dix2B(1234,10);B :-=dix2B(99999,10); 
S :=Somme (A,B) ;B2dix(S,10) 
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On obtient : 

[1,0,1,2,3,3],101233 

On vérifie : 1234 + 99999 = 101233 

Remarque Avec Xcas les polynômes peuvent être représenté sous forme symbol- 
ique (par exemple 11%? + 22x + 33) ou sous la forme de la liste de ses coefficients 
que l’on écrit par ordre décroissant (par exemple [11,22,33]). On peut passer d’une 
représentation à l’autre avec : symb2poly ou poly2symb. 

Par exemple : 

symb2poly(11xx”2+22+xx+33)=[11,22,33] et 

normal (poly2symb([11,22,33]))=11xx"2+22x+x+33. 

Si les listes À et B représentent les coefficients de 2 polynômes (coefficients que 
l’on écrit par ordre décroissant), alors pour faire la somme de ces 2 polynômes, 
l’algorithme est le même à condition de ne pas gérer les retenues (par exemple la 
liste [11,22,33] représente le polynôme 11x? + 22x + 33). 

On tape : 


//A et B sont les listes des coefficients de 
//2 polynomes (ordre decroissant) 
//on renvoie la liste A+B 
Sommep (A,B) :={ 
local j,sa,sb,s,c,Result; 
Result:=f[]; 
sa:-=dim(A)-1; 
sb:-=dim(B)-1; 
si sa>sb alors 
B:=concat ([0$ (sa-sb)],B);s:=sa; 
sinon 
A:=concat ([0$ (sb-sa)],A);s:=sb; 
fsi; 
ret:=0; 
pour j de s jusque 0 pas -1 faire 
G:=AlI)+BET; 
Result[jl:=c; 
fpour; 
retourne Result; 


}:; 


On tape par exemple : 

A :=[12,-3,35]; 
B :=[23,1,23,45]; 
S :=Sommep (A,B) 
On obtient : 
[23,13,20,80] 


9.3 Produit d’un entier par un entier a <10 


On va écrire la fonction produit d’un entier n par un entier a < 10 comme on le 
fait à la main. Pour cela, on va prendre comme paramètre de la fonction Prodan 
le nombre a et la liste N constituée par les nombres servant à l’écriture de l’entier 
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n (on pourra transformer le nombres n en cette liste avec dix2B et obtenir N). 
On tape : 


//on entre un entier a<10 et une liste N de nombres 
//compris entre 0 et 9 
//on renvoie une liste de nombres compris entre 0 et 9 
Prodan(a,N):={ 

local j,sn,c,ret,Result; 

si a>=10 alors retourne "erreur"; fsi; 

Result:=[]; 

sn:=dim(N)-1; 

ret:=0; 

pour j de sn jusque 0 pas -1 faire 

c:=axN[j]+tret; 


Result[jl:-irem(c,10); 
ret:=iquo(c,10); 
fpour; 


si ret>0 alors Result:=-prepend(Result,ret); fsi; 
retourne Result; 


On tape par exemple : 

a :=4;N :=dix2B(199999,10); 

Pa :-=Prodan(a,N) ;B2dix(Pa,10) 

On obtient : 

[7,9,9,9,9,6], 799996 

On vérifie : 4 x 199999 — 799996 

Remarque Si la liste P représente les coefficients d’un polynôme (coefficients que 
l’on écrit par ordre décroissant), alors pour faire la produit de ce polynôme par un 
réel a, l’algorithme est le même à condition de ne pas gérer les retenues. 

On tape : 


//on entre un entier a et la liste P des coefficients d’un 


//polynome (ordre decroissant) 
//on renvoie la liste axP 
Prodanpi(a,P):={ 
local j,sn,c,Result; 
Result:=[]; 
sn:=dim(P)-1; 
pour j de sn jusque 0 pas -1 faire 
c:=axPl[;j]; 
Result[jl:=c; 
fpour; 
retourne Result; 


}:; 


On tape par exemple : 
a :=11; 
P =123;1:23,451,; 
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Pa :-=Prodanpi{(a,P) 
On obtient : 
[253,11,253,495] 


9.4 Produit de 2 entiers 


On va écrire la fonction produit de deux entiers comme on le fait à la main. 
Pour cela on va prendre comme paramètre de la fonction Produit les listes À et 
B constituées par les nombres servant à l’écriture de ces deux entiers (on pourra 
transformer chacun des 2 nombres a et b en ces listes avec dix2B et obtenir À et 
B). On se servira des deux fonctions précédentes. 

On tape : 


//on entre des listes de nombres compris entre 0 et 9 
//on renvoie une liste de nombres compris entre 0 et 9 
Produit (A,B) :={ 
local j,sa,sb,C,Result; 
sb:=dim(B)-1; 
Result:=[0$sbl]l; 
pour j de sb jusque 0 pas -1 faire 
C:=Prodan(Bfj],A); 
C:=concat (C, [0$(sb-;j)]); 
Result :=Somme (Result, C); 
fpour; 
retourne Result; 
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On tape par exemple : 

A :=dix2B(1234,10);B :-=-dix2B(9999,10) ; 

P :-=Produit (A,B) ;B2dix(P,10) 

On obtient : 

[1,2,3,3,8,7,6,6],12338766 

On vérifie : 1234 x 9999 = 12338766 

Remarque Si les listes À et B représentent les coefficients de 2 polynômes (co- 
efficients que l’on écrit par ordre décroissant), alors pour faire le produit de ces 2 
polynômes, l’algorithme est le même à condition de ne pas gérer les retenues. 

On tape : 


//A et B sont les listes des coefficients de 
//2 polynomes (ordre decroissant) 
//on renvoie la liste qui represente le produit des 2 polynomes 
Produitp(A,B):={ 
local j,sa,sb,C,Result; 
sb:=dim(B)-1; 
Result :=[0$sb]; 
pour j de sb jusque 0 pas -1 faire 
C:=Prodanp(B[j]l,A); 
C:=concat (C, [0$(sb-j)]); 
Result :=Sommep (Result,C); 
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fpour; 
retourne Result; 
}:; 
On tape par exemple : 
A :=[12,-3,35]; 
B :=[23,1,23,45]; 
P :-=Produitp(A,B) 
On obtient : 
[276,-57,1078,506,670,1575] 
On vérifie et on tape : 
normal (poly2symb (A) xpoly2symb(B)) 
On obtient : 
276xx 5-57xx 4+1078xx 3+506xx 2+670xx+1575 


10 Jeu de recherche un nombre 


10.1 Le jeu 
10.2 C’est vous qui jouez 


L'ordinateur choisi au hasard un nombre n entre 1 et 1000 et vous devez le 
trouver. Chaque fois que vous proposez une réponse rep. l’ordinateur vous dit si 
le nombre est plus grand ou plus petit que votre réponse ou bien il vous dit "Bravo" 
en vous donnant votre nombre d’essais. 

On tape : 


je _trouve():={ 
local n,rep, j; 
n:=rand(1000)+1; 
saisir("trouvez un nombre entre 1 et 1000",repj); 
j:=1; 
tantque n!=rep faire 
j:=3j+1; 
si rep<n alors 
saisir("c'est plus grand que "+rep+", essai "+j,rep); 
sinon 


si rep>n alors 
saisir("c'est plus petit que "+rep+'", essai "+j,rep); 
fsi; 
fsi; 
ftantque; 
retourne "Bravo "+n+" trouve en "+j+" fois"; 
}:; 


Pour jouer on tape : je_trouve () 


10.3 C’est l’ordinateur qui joue 


Vous devez choisir un nombre n : on le fait avec choisir nombre. Quelque 
soit le nombre n, dans le programme dichot l’ordinateur choisit le milieu de 
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l'intervalle dans lequel se trouve ce nombre. Si l’ordinateur sait que le nombre 
choisi est entre m et M, l’ordinateur choisit floor ( (m+M) /2) caril doit proposer 
un nombre entier et la réponse proposée par l’ordinateur est stockée dans rep. 
Vous devez répondre : 

— "+" si votre nombre est plus grand que le nombre proposé, 
"-"°s1 votre nombre est plus petit que le nombre proposé, 
"=" si votre nombre est égal au nombre proposé. 
Votre réponse est mise dans la variable q : Si vous répondez n’mporte quoi on vous 
redemande une entée correcte : c’est ce que fait repeter. La variable : compte 
le nombre d’essai . On tape : 


choisir _nombre():={ 
local n; 
saisir("Choisissez un nombre entre 1 et 1000",n); 
afficher("l'ordinateur doit trouver ce nombre en vous interrogeant"); 
retourne n; 
}:; 
dichot () :={ 

local j,m,M,rep,q,n; 

m:=1; M:=1000; 

pour j de 1 jusque 11 faire 

rep:=floor ((m+M)/2); 


repeter 

saisir _ chaine "essai "++ ": "+rep+" est-ce +, —- ou =",q; 
jJusqua qg=="+" ou q=="-" où q=="="; 
si q=="=" alors 


afficher(rep+" trouve en "+j+" essais"); retourne rep; 
fsi; 
si q=="+" alors 
m:=rep+l; 
sinon 
M:=rep-l; 
fsi; 
si m—=M alors 
afficher(m+t" trouve en "++" essais"); retourne m; 
fsi; 
fpour; 
}:; 
jeu_ordi():={ 
local n,rep; 
n:=choisir nombre); 
rep:=dichot(); 
afficher("j’avais bien choisi :"+n); 
retourne rep; 


fie 


On tape : 
jeu_ordi () 
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10.4 Solution approchée de f(x) = 0 sur [a; b] 


On suppose que f est continue et que f(a) x f(b) < 0. Donc f s’annule sur 
[a; b]. On va procéder par dichotomie : on partage l’intervalle [a; b] en deux : 


La, b] = [a;m = (a +b)/2] U [m = (a +b)/2; bd] 


et on regarde si f(a) x f(m) < 0 si c’est le cas on cherche la solution dans [a; m] 
sinon on cherche la solution dans [m; b]. On recommence autant de fois qu’il faut 
pour avoir la précision eps désirée. 

On tape : 


Zero(f,a,b,eps):=1{ 
local m; 
a:=evalf(a); b:-=evalf(b); 
si f(a)==0 alors retourne a fsi; 
si f(b)==0 alors retourne b fsi; 
si f(a)xf(b)>0 alors retourne "erreur" fsi; 
tantque (abs(a-b)>eps) faire 
m:=(a+b)/2; 
si (f(a)x£f(m)<0) alors 
b:=m; 
sinon 
a:=m; 
fsi; 
ftantque; 
si a<b alors retourne aàa,b; sinon retourne b,a; fsi; 


PS 


On tape : 
Zero(x->x"2-2,0,2,1e-10) 

On obtient : 
1.41421356233,1.41421356238 


11 Les carrés emboités 


On veut tracer n carrés : un carré direct de côté AB = d, puis les carrés directs 
de côtés AB; lorsque les B; définissent une subdivision de AB en n parties égales 
(AB: =d/n=BiBs=.:=8B, 1B). 

11.1 Les carrés emboités et la tortue 
11.1.1 Les carrés emboités (itératif) 


On tape dans un éditeur de programmes : 


Carresl(d,n) :={ 
local j,k; 
k:=d/n; 
pour j de 1 jusque n faire 
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repete(4,avance d,tourne gauche 90); 
d:=d-Kk; 
fpour; 
}:; 


11.12 Les carrés emboités (récursif) 


La figure Carrerl(d,n) estcomposée du grand carréetdeCarrerl(da-d/n,n) 
On tape : 


Carrerli(d,n):=f{ 
si n<0 alors retourne 1;fsi; 
repete(4,avance d,tourne_ gauche 90); 
d:=d-d/n; 
Carrerl(d,n-1); 


r 

On ouvre un écran de géométrie tortue et on tape avec Alt +d : 
efface; 

Carres1(100,6); 

pas_de_ cote -120; 

Carrer1(100,6); 


11.2 Les carrés emboités en géométrie 
11.2.1 Les carrés emboités (itératif) 
On tape : 


carres(A,B,n):={ 

local j,L,k; 

L:=NULL; 

k:=(B-A)/n; 

pour j de 1 jusque n faire 
B:=translation(k+xj,A); 
L:=L,carre(A,B); 

fpour; 

retourne L; 


}:; 

Puis on tape dans une ligne d’entrée : 
carres (point (0),point (100),6) 
et on voit le dessin dans l’écran réponse. 


11.2.2 Les carrés emboités (récursif) 


On tape dans un éditeur de programmes : 


carrer(A,B,n):={ 
si n<0 alors retourne 1;fsi; 
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carre (A,B); 
B:=translation((B-A)x(n-1)/n,A); 
retourne carrer(A,B,n-1); 


}:; 


Puis on tape : 
DispG() ; carrer(point (0),point (100),5) 
et on voit le dessin dans l’écran DispG car c’est seulement ce qui est retourné par 
la fonction qui est dessiné dans la réponse, alors que toutes les objets graphiques 
intermédiaires sont affichés dans l’écran DispcG. 

Ou bien on tape dans un éditeur de programmes : 


carrerL(A,B,n,L):={ 
si n<0 alors retourne L;fsi; 
L:=[L,carre(A,B) |; 
B:=translation((B-A)x(n-1)/n,A); 
retourne carrerL(A,B,n-1,L); 


}:; 


Puis on tape dans une ligne d’entrée : 
carrerL(point (0),point (100),6,[]) 
et on voit le dessin dans la réponse. 

Remarque pour les utilisateurs avancés 
On peut améliorer le programme précédent en utilisant une seule liste que l’on 
modifie en place (avec l’opérateur =<) afin de ne pas recopier la liste L à chaque 
affectation par := (un peu comme avec un passage d’arguments par référence). 
Attention, toutefois, cela nécessite de faire une copie de la liste vide initiale par 
copy sinon c’est la liste du programme lui-même qui sera modifiée et ne sera 
donc plus initialisée à une liste vide. 


carreenplace(A,B,n,L):={ 
si n<0 alors retourne L;fsi; 
L[dim(L)] =< carre(A,Bj); 
B:=translation((B-A)x(n-1)/n,A); 
retourne carreenplace(A,B,n-1,L); 
}:; 


appel _carre(A,B,n):={ 
local L; 
L:=copy([]l); 
carreenplace(A,B,n,Lj); 
retourne L; 


}:; 


12 D’autres carrés emboités 


On veut tracer n carrés : on trace un carré direct de côté AB = d, puis le carré 
de sommets les milieux des côtés du carré dessiné précédemment... 
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12.1 Avec la tortue 


12.1.1 Les carrés (itératif) 


On tape : 


Carres21(d,n):={ 


}:; 


local j; 
pour j de 1 jusque n faire 


repete 4,avance d,tourne gauche 


avance d/2; 
tourne _ gauche 45; 
d:=d/sqrt (2.); 

fpour; 

pour j de 1 jusque n faire 
d:=dxsqrt(2.); 
tourne droite 45; 
avance -d/2; 

fpour; 


LA 


90; 


On tape si on ne veut pas repasser par les mêmes traits : 


Carres22 (d,n):={ 
local c,j,k; 
Crs=d; 
j:=1; 


}:; 


tantque j!=n+1 faire 
avance c/2; 
tourne _ gauche 45; 
ci=c/sqrt(2.); 
j:=3+1; 
ftantque; 
tantque j!=1 faire 
tourne droite 45; 
c:=cxsqrt(2.); 
avance c/2; 
k:=1; 
tantque k!=-4 faire 
tourne_gauche 90; 
avance c; 
k:=k+1; 
ftantque; 
tourne _gauche 90; 
Jde 
ftantque; 


LA 


12.1.2 Les carrés (récursif) 


On tape : 


27 


Carrer21(d,n):={ 
si n<=0 alors retourne 1; fsi; 
repete 4,avance d,tourne_ gauche 90; 
avance d/2; 
tourne_gauche 45; 
Carrer21(d/sqrt(2.),n-1); 
tourne droite 45; 
recule d/2; 

}:; 


On tape si on ne veut pas repasser par les mêmes traits : 


Carrer22(d,n):={ 
si n>0 alors 
avance d/2; 
tourne _ gauche 45; 
Carrer22 (d/sqrt(2.),n-1); 
tourne droite 45; 
avance d/2; 
tourne_gauche 90; 
repete 3,avance d,tourne gauche 90; 
fsi; 


hic 


On ouvre un écran de géométrie tortue et on tape avec Alt +d : 
efface; 

Carres21(100,6); 

saute -120; 

Carres22(100,6); 

pas_de_ cote -120; 

Carrer21(100,6) ; saute 120; 
Carrer22(100,6) ; 


12.2 Les carrés en géométrie 
12.2.1 Les carrés (itératif) 
On tape : 


carres2(A,B,n):={ 

local j,L,C; 

L:=NULL; 

pour j de 1 jusque n faire 
L:=L,carre(A,B,C); 
A:=point ((A+B)/2); 
B:=point ((C+B)/2); 

fpour; 

retourne L; 


}:; 
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Puis on tape dans une ligne de commande : 
carres2 (point (0) ,point (100),6) 
et on voit le dessin dans l’écran réponse. 


12.2.2 Les carrés (récursif) 
On tape : 


carrer2(A,B,n):={ 
local C; 
si n==1 alors retourne carre(A,B) 
carre(A,B,C); 
A:=point ((A+B)/2); 
B:=point ((C+B)/2); 
retourne carrer2(A,B,n-1); 


fire 


Puis on tape : 
carrer2 (point (0),point (100),6) 
et on voit le dessin dans l’écran Di spG. 

Ou bien on tape : 


carrerL2(A,B,n,L):={ 

local C; 

si n==0 alors retourne L fsi; 

L:=[L,carre(A,B,cC)l; 

A:=point ((A+B)/2); 

B:=point ((C+B)/2); 

retourne carrerL2(A,B,n-1,L); 
}:; 


On tape alors dans une ligne de commande : 
carrerL2 (point (0),point (100),6,[]) 
et on voit le dessin dans l’écran réponse. 
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fsi; 


